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1.  Introduction;  At  present  there  Is  widespread  Interest  In  the  stable  adaptive 


control  of  unknown  linear  time -Invariant  plants  using  Input-output  data.  Schemes 
have  been  suggested  for  both  direct  [1-3]  and  Indirect  [4,5]  control  of  continuous 
as  well  as  discrete  [6,7]  systems  and  the  equivalence  of  the  two  schemes  In  some 
cases  has  also  been  demonstrated  [4,5].  Probably  the  single  most  Important  prob- 
lem to  arise  In  the  course  of  these  Investigations  concerns  the  proof  of  stability 
of  the  overall  adaptive  control  loop. 

Monopoll  [1]  proposed  a scheme  for  continuous  systems  Involving  an  auxiliary 
signal  fed  Into  the  reference  model  and  a corresponding  augmented  error  between 
model  and  plant  outputs.  Narendra  and  Valavanl  [2],  using  p< sltlve  real  operators, 
suggested  a similar  approach  and  clarified  the  resulting  stability  problem  when 
the  relative  degree  of  the  plant  Is  greater  than  or  equal  to  three.  They  offered 
a conjecture  that  the  adaptive  loop  would  also  be  stable  for  the  general  case. 

Feuer  and  Morse  [3]  proposed  a stable  solution  to  the  adaptive  control  problem  but 
the  resulting  controller  Is  much  too  complex  for  use  In  practical  applications. 
Thus,  the  search  has  continued  for  a controller  with  a simple  structure  which  will 
assure  the  asymptotic  stability  In  the  large  of  the  adaptive  loop.  The  results 
presented  In  this  paper  demonstrate  the  desired  stability  behavior  for  discrete 
versions  of  the  simple  controllers  suggested  In  [1]  and  [2].  The  control  problem 
^^^^he  continuous  case  however  remains  unresolved. 

This  paper  examines  the  discrete  version  of  the  problem^ considered  In  [2] 


^recapitulating  the  basic  philosophy  as  well  as  the  specific  technique  used  for 
the  design  of  the  adaptive  controllei;^l!^that  paper':  Hence^he  first  few  sections 
of  this  paper  have  been  considerably  condensed  and  the  interested  reader  Is  referred 
to  the  earlier  work  for  all  details.  The  principal  contribution  made  here  is  the 
verification  of  the  conjecture  mads  [ij^regarding  the  stability  of  the  adaptive 
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loop,  for  the  discrete  problem.  Accordingly  most  of  the  paper  Is  devoted  to  the 
proof  of  stability.  These  results  together  with  some  recent  work  done  on  the 
stability  of  differential  equations  with  unbounded  coef flclents-+8^may  prove 
helpful  In  resolving  the  stability  problem  of  continuous  adaptive  systems  as  well.^^ 
2.  Statement  of  the  Problem; 

A single-input  single-output  discrete  linear  time-invariant  plant  P Is  described 

by  the  Input-output  pair  {u(k),y  (k)}  and  can  be  represented  by  the  transfer  function 

P 


Ziz) 

p 


(1) 


where  W (z)  Is  proper  with  R (z)  a monlc  polynomial  of  degree  n,  Z (z)  a monlc 
P P P 

stable  polynomial  of  degree  m < n and  k a constant  gain  parameter.  The  Integer 

P 

n - m Is  called  the  relative  degree  of  the  plant.  We  assume  that  only  m,n  and  the 

sign  of  k are  known  while  the  coefficients  of  Z and  R are  unknown. 

P P '^P 

A reference  model  M whose  output  y^Ck)  represents  the  behavior  desired  from 
the  plant  when  augmented  by  a suitable  controller  can  be  represented  by  the  transfer 
function 


(2) 


where  monlc  stable  polynomials  of  degrees  n and  r < m respec- 

tively and  k^  Is  a constant.  Hence  the  relative  degree  of  the  model  Is  assumed  to 
be  greater  than  or  equal  to  that  of  the  plant.  The  Input  to  the  model  Is  a speci- 
fied uniformly  bounded  reference  Input  r(k) . 

The  adaptive  control  problem  Is  to  determine  a suitable  control 


such  that 


e.Ck)  « y„(fc)  - ■*■0  as  k 


with  all  zeros  Inside  the  unit  circle. 
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For  the  sake  of  simplicity  we  shall  assume  that  r - m.  As  In  the  continuous 
case  the  solution  to  the  above  problem  may  be  divided  Into  two  parts.  The  first 
part  which  Is  algebraic  In  nature  addresses  Itself  to  the  realizability  of  a 
suitable  controller  structure.  It  can  be  shown  exactly  as  In  the  continuous  case 
[2]  that  a controller  can  be  found  which  can  achieve  (3)  with  a fixed  set  of 
parameters.  In  the  following  section  the  equations  describing  the  controller  are 
merely  stated.  The  second  part  Is  analytic  In  nature  and  deals  with  the  stability 
of  the  adaptive  error  equations.  Again,  It  Is  found  that  when  m ••  n (or  m-n-1  In 
the  continuous  case)  the  adaptive  equations  can  be  shown  relatively  easily  to  be 
asymptotically  stable.  Hence  our  main  Interest  Is  In  the  case  m < n - 1 and 
auxiliary  Inputs  have  to  be  fed  Into  the  reference  model.  The  statement  of  the 
stability  problem  and  the  proof  of  stability  are  considered  In  detail  In  section  5. 
3.  Structure  of  the  Adaptive  Controller; 

As  in  12]  two  different  structures  are  needed  for  the  discrete  adaptive  con- 
trol problem  corresponding  to  the  two  cases  m ■ n and  m < n - 1.  When  m « n and 

* 

the  model  transfer  function  Is  assumed  to  be  positive  real  the  simple  structure 
shown  In  Figure  (1)  can  be  used.  For  the  case  when  m < n - 1,  as  described  In  [2], 
an  auxiliary  signal  has  to  be  fed  into  the  model  and  the  corresponding  structure 
Is  shown  In  Figure  (2) . A brief  description  of  the  controller  structure  for  the  two 
cases  Is  given  below: 

Case  (1)  (m*n) t The  controller  consists  of  (2n+l)  adjustable  parameters  which  are 
the  elements  of  a parameter  vector  e(k)  defined  by 

e^(k)  - [kQ(k),c^(k),...c^(k),d^(k),d2(k),...d^(k)]. 


Two  Identical  auxiliary  signal  generators  of  dimension  'n*  having  state  variables 
v^^\k)  and  v^^^Ck)  and  Inputs  u(k)  and  y^Ck)  respectively  as  shown  In  Figure  (1) 


* 

There  Is  no  loss  of  generality  here  since  by  prefllterlng  the  model  can  be  made 
positive  real. 


r 
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, (1)/unT  (2).,  .T, 

w (k)  - [v  (k)  ,v  (k)  1 

C(k)  - L"^(B)w(k) 


Figure  2 
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form  part  of  the  controller.  If  a vector  w(k)  is  defined  as 


- [r(k).v^^^(k)'^,v^^\k)’'l 


the  signal  fedback  Into  the  plant  may  be  represented  by 


e^(k)a)(k)  - aii;(k)’^rw(k)ej^(k) 


The  first  term  In  (4)  represents  a linear  combination  of  the  elements  of  (o(k)  and 
corresponds  to  the  feedback  signal  In  the  continuous  case.  The  second  term  which 
depends  on  the  output  error  e^(k)  Is  found  essential  In  the  discrete  case  to 
establish  the  stability  of  the  error  equations  as  described  In  [9] . If  u(k)  Is 
bounded,  as  the  adaptation  proceeds,  this  term  Is  seen  to  tend  to  zero  with  e^(k). 

—It 

Following  the  results  In  I2J  It  can  be  shown  that  a constant  vector  6 exists 
such  that  %rhen  9(k}  ^ 9 the  transfer  function  of  the  plant  together  with  the 
controller  matches  that  of  the  reference  model. 

The  adaptive  control  problem  Is  to  determine  the  law  for  updating  e(k)  such 

that  0(k)  '►0  as  k “ while  maintaining  overall  system  stability. 

Case  Cil)  Gn^n-1) ; With  no  loss  of  generality  we  can  assume  that  L(z)  a rational 

function  In  z (with  L ^(z)  a strictly  proper  minimum  phase  function)  exists  such 

that  W (z)L(z)  Is  strictly  positive  real.  However,  since  L(z)  Is  not  physically 

realizable  the  same  modification  as  that  suggested  In  the  continuous  case  has  to 

be  used  here  as  well.  As  shown  below,  this  Involves  feeding  back  signals  Into  both 

plant  and  model  such  that  the  error  equations  have  the  same  form  as  In  case  (1) 

(ref.  Section  4).  Figure  (3)lndlcates  the  structure  of  the  controller  when  the  plant 

gain  k Is  known;  for  simplicity  It  Is  assumed  that  k ■ k^  - 1.  Since  In  this 
P p M 

T 

case  only  2n  parameters  have  to  be  adjusted  we  define  6 (k)  > [c^(k) •C2(k} , . . . 

. . ,c^(k)  ,dj^(k)  , . . . ,d^(k)  ] and  ui^CW  - tv^^\k)^,v^^^  (k)^] . 

It  is  obvious  that  a rational  function  L(z)  with  denominator  polynomial  Z^(z) 
exists  such  that  W^(z)L(z)  is  strictly  positive  real. 


The  signal  fed  back  Into  the  plant  Is 


u (k)  - 9^(1)  (k) 

P 

while  the  signal  fed  back  Into  the  model  Is 

Ujj(k)  - L(z)  [{L”^(z)0’^(k)-0’^(k)L"^(2)}a)(k)  + ac’^(k)rc(k)e^(k)  ] 

(where  L ^(z)u)(k)  » C(k)) 

so  that  the  resulting  error  equations  have  the  form  required  to  generate  stable 

adaptive  laws  as  described  In  Section  4. 

While  L(z)  Is  not  physically  realizable,  W„(z)L(z)  can  be  realized  and  hence 

n 

the  overall  system  Is  as  shown  In  the  figure  (3). 

I y desired  output 
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4.  The  Error  Equations; 

Let  the  parameter  vector  6(k)  be  expressed  as 

e(k)  - e*  + ♦(k) 


where  4'(k)  represents  the  parameter  error  vector  at  time  k. 
Case  (1)  (m«n)  ; The  output  error  e^^Ck)  may  be  expressed  as 

— T — — T — 

(k)a)(k)  - oo)  (k)  r(i)(k)ej^(k) } 


where  Is  a strictly  positive  real  operator.  From  the  recent  results  In  [9] 

(also  given  In  detail  In  the  next  section)  It  Is  seen  that  If  ^^Ck)  Is  updated 
according  to  the  law 


Kk+1)  - ?(k)  - rej^(k)a!(k) 


"*■  0 whether  or  not  w(k)  Is  bounded.  Since  In  this  case  the  (desired)  output 
of  the  model  Is  uniformly  bounded,  the  plant  output  will  also  be  uniformly  bounded 
and  approach  the  desired  output  asymptotically. 

Case  (11)  Cm<n-1) ; The  output  error  ej^(k)  In  this  case  (see  [2])  satisfies  the 
error  equation 


ei(k)  - [Wj^(z)L(z)]{/(k)i;(k)  - oi:’^(k)rc(k)ej^(k) } 

where  W (z)L(z)  Is  a strictly  positive  real  transfer  function  and  as  defined  earlier 
M 

L‘^(z)u)(k)  - C(k). 

By  the  same  arguments  as  in  case  (1)  it  follows  easily  that  If  the  adaptive  law 


is  used,  then 


♦ CW-1)  - - re^(k)ak) 


ej(k)  -►0  as  k 


However,  It  no  longer  follows  that  the  plant  output  will  be  bounded  since  the  model 

* 

* This  applies  to  case  (i).  For  case  (il)  6(k)  " & + ^(k). 
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output  (which  is  due  to  both  reference  and  the  auxiliary  inputs)  may  be  unbounded. 
Hence,  to  prove  the  global  asymptotic  stability  of  the  adaptive  system  it  is 
necessary  to  show  that  neither  the  plant  output  nor  the  model  output  can  be  un- 
bounded - in  other  words  verify  for  the  discrete  case  the  conjecture  made  in  [2] 
for  continuous  systems.  The  rest  of  this  paper  is  devoted  entirely  to  this  problem. 
5.  Verification  of  the  Conjecture  by  Narendra  and  Valavani  for  Discrete  Systems; 
a)  Description  of  the  Error  Model; 

The  error  model  whose  stability  is  to  be  analyzed  is  a complex  vector  nonlinear 
difference  equation.  For  convenience  of  analysis  we  shall  consider  it  in  three 
separate  parts  which  correspond  to  the  three  parts  of  the  system  shown  in  Figure(4) . 


Part  I Part  II  Part  III 

The  Plant  Feedback  Loop  The  Prefllter  The  Error  Model 

Figure  4 

Part  I - The  Plant  Feedback  Loop;  The  plant  together  with  the  controller  can  be 
described  by  the  vector  difference  equation 


x(Irt-l)  - A xCk)  + b„H^(k)uCk)  + r(k)] 

B B 

u(k)  - C_x(k) 

B 


(5) 
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a)(k)  represents  the  output  vector  of  Interest,  x(k)  is  the  state  vector  of  the 

plant  together  with  the  controller  and  when  <ti(k),  the  parameter  error  vector  is 

identically  zero,  the  plant  and  model  transfer  functions  match  exactly.  r(k),  the 

reference  input,  is  uniformly  bounded  and  the  matrices  A and  C , and  vector  b 

are  of  appropriate  dimensions.  As  described  in  the  previous  section  A^  is  a 

(3nx3n)  stable  matrix  and  b is  a (3nxl)  vector  and  C is  a (2nx3n)  matrix. 

m m 

Part  II  - The  Prefilter; 

The  second  part  of  the  system  shown  in  Flgure( 4) consists  of  a diagonal  trans- 
fer matrix  all  of  whose  elements  are  the  same  and  equal  to  W (z)  » L ^(z)  or  ' 

Li 

C^Ck)  - L‘^(z)u)^(k)  (6) 

L ^(z)  is  assumed  to  be  an  asymptotically  stable  system  of  relative  degree  n - m as 
described  earlier  (number  of  poles  - number  of  zeros  » n - m)and  of  minimum  phase. 

Part  III  - The  Error  Model  [91; 

The  third  part  is  the  model  of  the  error  equations  described  in  the  previous 

section  and  consists  of  a strictly  positive  real  transfer  function  in  the  feed- 

T T 

forward  path  and  feedforward  and  feedback  gains  ^ (k)  and  ac  (k)rc(l()  respectively 
as  shown  in  Figure(4).  It  can  also  be  represented  by  a 3n  order  difference 
equation 


e(lcfl)  “ A e(k)  + bv(k) 

ID 

ej^(k)  - c^e(k)  + dv(k)  (7) 

v(k)  - ♦’'(k)i;(k)  - ac’'^(k)rc(k)e^(k) 

1 T 

a > y,  r - r > 0 

T -1 

where  d -f  c (zI-A^)  b is  strictly  positive  real. 

The  pereaetcr  error  vector  ^(k)(and  hence  the  parameter  vector  e(k))  is  ad- 


justed according  to  the  law 

♦ (k+1)  - ♦(«  -re^(k)C(k) 
(or  e(lcfl)  - e(k)  - rej^(k)C(k)) 


(8) 
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The  four  sets  of  difference  equations  (5), (6), (7)  and  (8)  completely  determine  ] 

« 

the  error  model  of  the  overall  discrete  system.  The  signal  C(k)  and  output  error  . 

1 

e^(k)  determine  how  (^(k)  Is  updated  but  this,  In  turn,  determines  the  nature  of  <»(k)  I 

and  C(k) . . 

b)  Statement  of  the  Conjecture;  The  conjecture  in  [2]  when  applied  to  the  problem 
above  may  be  stated  as  follows: 

If  <^(k)  Is  adjusted  according  to  the  law  (8)  to  keep  the  output  error  bounded,  ‘ 

I 

the  output  of  the  plant  l.e.  a)(k)  will  also  be  bounded. 

Equivalently,  Che  conjecture  Implies  that  the  overall  nonlinear  system  described  j 

by  (5),C6),C7)  and  (8)  Is  stable  and  that  all  the  signals  are  uniformly  bounded.  ^ 

Since  the  stability  In  the  large  of  the  above  nonlinear  system  Is  Intractable,  I 

J 

we  shall  consider  the  three  linear  blocks  In  Flgure( 4) separately  to  simplify  l 

1 

analysis.  ^ 

c)  A Qualitative  Analysis;  , 

Fart  I of  Flgure(4) Is  a feedback  loop  with  a stable  time-invariant  forward 
path  and  a time-varying  gain  vector  i^(k)  In  the  feedback  path.  The  output  vector 
(o(k}  can  be  either  uniformly  bounded  or  unbounded.  In  the  former  case  ^(k),  the 
output  of  the  pref liter  Is  also  uniformly  bounded  and  the  behavior  of  part  III 

it 

and,  hence,  chat  of  the  entire  system  are  completely  known.  If  Che  Input  r(k)  j 

Is  sufficiently  rich,  the  vectors  <i>(k)  and  C(k)  will  be  uniformly  bounded  while 
^(k)  and  e^Ck)  tend  to  zero  as  k If,  however,  r(k)  Is  not  sufficiently  . 

rich,  ^(k)  , a constant  vector,  while  ej^(k)  0.  In  both  cases  the  signals  ^ 

In  the  system  are  bounded.  If,  however.  It  Is  assumed  that  u(k)  and,  hence  C(k)  t 

i 

are  unbounded,  Che  analysis  In  the  following  section  shows  that  we  are  led  to  a k 

contradiction.  Hence,  only  the  first  alternative  Is  possible  (l.e.  u)(k)  Is  | 

uniformly  bounded)  and  the  conjecture  Is  verified.  ( 

4 
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Before  proceeding  to  give  an  analytic  proof,  ve  present  here  a brief 
({ualltatlve  analysis  of  the  various  steps  Involved. 

In  section  d.  It  Is  first  shown  that  e^(k)  and  e^(k)C(k)  tend  to 

zero  asymptotically  whether  or  not  C(k)  Is  uniformly  bounded.  This,  In 

* 

turn.  Implies  that  ^(k)  tends  to  a constant  vector  ♦ . The  asymptotic  be- 
havior of  part  I can  be  described  In  terms  of  the  dominant  modes  of  the 

. ip 

unstable  system  x ■ [A^  + b^j^^  (k)C^]x  and  It  Is  shown  that  asymptotically 
b>(k)  and  C(k)  exhibit  similar  behavior  and  that  ^(k)  is  asymptotically 
orthogonal  to  the  dominant  mode  of  C(k)  and  hence  io(k).  This,  In  turn. 
Implies  that  the  dominant  mode  of  <tf(k)  Is  not  present  In  the  feedback  signal 
In  part  I as  k which  contradicts  the  assumption  that  (i>(k)  Is  unbounded. 

Hence,  the  overall  system  Is  such  that  when  there  Is  no  reference 
signal  “(k),  C(k)  and  e^^Ck)  0 as  k * and  ♦(k)  -*■  4'*  a constant  vector 
such  that  the  feedback  loop  Is  asymptotically  stable.  When  r(k)  is  suffi- 
ciently rich,  <^(k)  and  ^(k)  are  uniformly  bounded  but  4>(k)  and  e^(k)  0. 


The  Error  Model; 


From  the  discrete  version  of  the  iC.''.lman-Yacubovlch  lemma  If 

T -1  T 

d c (zI-A  ) b Is  strictly  positive  real  a matrix  P > P >0  and  a vector 

n 

q exist  such  that 

A^^PA  - P - -qq^  - eL;  A ^Pb  - c/2  + vq;  d - b’’Pb  - (9) 

■ ID  ID 


I 

i 

j 

i 


-11- 


T 

for  some  L ■ L >0  and  scalars  e,v  > 0. 

If  a Lyapunov  function  candidate  for  the  set  of  equations  (9)  Is  chosen  as: 

V[e(k),*(k))  - V(k)  - 2e(k)'^Pe(k)  + ♦’^(k)r"^^(k) 

It  was  shown  In  [ 9]  that  dv(e(k) ,^(k))  ■ dv(k)  V(kfl)  - V(k)  luiy  be  expressed  as: 

ii7(k)  - -2[e^(k)q  - w(k)]^  - 2ee’’(k)Le(k) 

+ (l-2o)c^(k)rc(k)eJ(k) 

<0  If  a > 

Hence  the  systcn  Is  stable  and  e(k)  and  ^(k)  are  bounded  If  e(0)  and  ^(0)  are  bounded. 
Furthemore,  since 

m 

I E AV(k)|  - |V(-)  - V(0)1  < - 
k-0 

m 

I 2(e^(k)q  - w(k)]^  + 2ee^(k)Le(k) 
k-0 

+ (2o-l)5*(k)rc(k)eJ(k)  < - 

we  conclude  that 

e(k)  -*•  0 J(k)rc(k)eJ(k)  -*•  0 (10) 

and  Sj^Ck)  ■►0  as  k • 

eheebsr  or  not  c(k)  i*  unlfondy  bounded. 

Since  ■ -rej^(k)c(k) • from  (10)  we  have  A^(k)  *^0  as  k '*’  **.  Hanee  S(k)  Is  a 

Cauchy  sequence  uhlch  belonfs  Co  a conpact  sat  and  hence 

^ . * 

S(k)  S a constant  vector 


(11) 
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Further  slnc6  e^^Ck)  ■ c e(k)  + dv(k) , the  Input  v(k)  to  the  strictly  positive 
real  transfer  function  also  tends  to  zero.  Hence  we  have 

♦’^(k)C(k)  - oJ(k)n:(k)ej^(k)  -►  0 as  k -►  » ( 

If  the  vector  C(k)  Increases  without  bound  there  exists  a subsequence  C(k)  such 
that  ||C(k)(|  Increases  monotonlcally  with  k and  ♦'*'(k)C(k)  - aC^(k) rcCWe^^Ck)  -►  0. 
Since  ((k)e^(k)  -*■  0 (eqn.  (10))  It  follows  that 

♦'^(i£)c(iS)  - ^i(^^)  ||c(^^)ll  ( 


where 


i\,  i\, 

li(k)  -*•0  as  k -►  *. 


We  now  consider  the  Implications  of  equations  (12)  and  (13)  for  different 


of  Inputs  C(k). 


Case  1 - ♦(k)  and  C(k)  Scalars; 


From  (13)  It  follows  that  If  |c(k)!  > ••  then  ^(k)  -*■  0 and  since  ^(k)  has  a 
limit  >(k)  0 as  k The  Importance  of  this  lies  In  the  fact  that  the  theoretical 

arguments  used  In  the  following  more  general  cases  are  In  essence  similar  to  that 
used  here. 

In  the  plant  feedback  loop  If  #(k)  -*■  0,  since  A Is  an  asymptotically  stable 

fli 

matrix,  x(k)  and  hence  w(k)  and  C(k)  are  bounded  which  contradicts  the  assumption 
made  that  c(k)  Is  unbounded. 

Case  2:  c(k)  - Cj^(k)Iij^  + o(l)]  Is  a constant  vector. 


We  consider  here  a special  but  Important  case  of  that  treated  In  Case  3. 

a constant  direction  and  c(k)  asymptotically  approaches  this  direction.  By  (11) 

s 

the  parameter  error  vector  ^(k)  converges  to  a constant  vector  ^ and  If 

m 

^|A^(k)||<  <•  , from  the  results  on  almost  constant  systems  [10]  we  note  that 
both  w(k)  and  ((k)  exhibit  this  type  at  behavior. 
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From  (13)  we  have 


and  hence 


Tf\^  f\, 

t « V > ■ V 


||C(k)|l 


•♦■0  as  k “ 


This,  In  turn.  Implies  (following  Che  same  arguments  as  In  Case  3)  that  the ^ 

T 

feedback  signal  ^ (k)w(k)  does  not  contain  asymptotically  the  dominant  component 

of  w(k) . Hence  x(k),  w(k)  and  c(k)  cannot  be  unbounded. 

2n 

Case  3:  c(k)  ■ £ C.(k)t.  (1.  are  2n  Independent  constant  vectors). 

1-1  ^ ^ ^ 

In  this  case,  we  consider  a general  unbounded  Input  vector  ((k)  whose 

"dominant  components"  lie  In  an  m^-dlmenslonal  subspace. 

Let  C,(k),(  (k),...  c (k)  be  unbounded 
1 m2 

while  ( (k),...,  C.  (k)  are  uniformly  bounded. 

m2'*'-*-  *n 

Further,  among  the  m2  unbounded  components  we  assume  Chat  m^(^m2)  are 
dominant  components  such  Chat 

supU  (v)|  -v-  sup||C(v)j|  l-l,2,...,m. 

k^v  Ittv 

I <^j^(k)|  - o[sup  ||C(v)  II  ] 1 - mj^+l,...,m2 

[Refer  to  Appendix  for  definition  of  o(*)  and  'v]. 

We  further  assume  (with  no  loss  of  generality)  that  C^(k)  are  "asymptotically 
Independent"  as  defined  In  the  appendix. 

e A 

Our  aim  here  Is  to  show  that  ^(k)  ^ such  that  ^ Is  orthogonal  Co  the 

subspace  generated  by  I,,!. and  hence  to  the  dominant  components  of  c(k). 

since  Che  dominant  components  of  w(k)  lie  In  the  same  subspace  If  W (z)  Is  a 

L 

* 

minimum  phase  transfer  function  ^ Is  also  orthogonal  to  these  components  and 
this.  In  turn,  leads  to  a contradiction  If  w(k)  Is  unbounded.  [Mote  that  If 
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~ ■ 1 the  present  case  reduces  to  Case  2] . 

By  (12) 

n m 

(k)A Jc  (k)  - olsup  ||c(v)||  ] (15) 

1-1  ^ ^ k^v  " 

Since  C^(k)  are  asymptotically  linearly  Independent  It  follows  that 

T 

11m  /(k)£  - (^*  1 - 0 (16) 

lr+»  ^ ^ 

2n 

If  w(k)  ■ E w (k)l  It  follows  from  (6)  that 
1-1  ^ 

C^(k)  - Wj^(z)w^(k)  l-l,2,...2n. 


Since  W (z)  Is  a strictly  proper  minimum  phase  transfer  function  It  follows  that 
L 

w^(k)  are  unbounded  for  1 - l,2,...,m2  and  uniformly  bounded  for  1 - m2^^,...,2n 

and  that  w. (k) , . . . ,w  (k)  are  the  dominant  components  of  w(k) . 

1 

The  feedback  loop  In  I Is  described  by  the  difference  equation 


x(k+l)  - [A  + b (|i^(k)C  ]x(k)  + b r(k) 

Q m m m 

- A x(k)  + b (()^(k)w(k)  + b r(k)  (17) 

n ID  ID 


since  r(k)  Is  uniformly  bounded  and  A Is  an  asymptotically  stable  matrix  It  follows 

m 

that  If  w(k)  Is  unbounded 


|x(k)|  < c sup  U^(v)w(v)|  + c,  ®1»®9  ^ 

^ k^v  ^12 

By  (16)|^^(k)w(k)|  - o[sup  ||  w(v)  ||  ] 

k^v ' 

and  hence 

|x(k)|  - o[sup||C  x(v)  II  1 
kiv  ■ 

^ o[sup||x(v)  II  1 (19) 

kiv'  ^ 

which  is  a contradiction  If  x(k)  Is  unbounded. 


i 
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Hence  x(k) , w(k)  and  C(k)  are  uniformly  bounded  and  the  adaptive  control  system 
Is  stable.  Further,  the  auxiliary  Input  to  the  model  Is 

u^(k)  - L(z){[L‘^(z)*(k)  - (^(k)L‘^(z)]w(k)  + ae^(k)c^(k)rc(k)} 

Since  ^(k)  and  e.  (k)  ■+•  0 as  k -►  »,  u (k)  -►  0 as  k •+■  “>  and  hence  the  model 

1 M 

output  approaches  the  desired  output  asymptotically.  Further  since  e^(k)  0 

the  plant  output  also  asymptotically  approaches  the  desired  output. 

6.  Conclusion;  The  paper  presents  a proof  of  the  stability  of  the  adaptive 
control  system  suggested  by  Narendra  and  Valavanl  [2]  for  the  discrete  case.  The 
same  proof  could  also  be  used  to  show  that  the  discrete  controller  suggested  by 
lonescu  and  Monopoll  [6]  Is  also  stable.  Hence  discrete  adaptive  systems  which 
use  an  augmented  error  signal  are  now  practically  feasible. 

The  adaptive  control  problem  In  the  continuous  case  using  an  augmented  error 


signal  still  remains  unresolved.  Recent  results  given  In  [8]  used  In  conjunction 
with  the  arguments  used  In  this  paper  may  prove  effective  In  Its  resolution. 
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Appendlx 

Definition  l!  Let  {x(k)}  and  {y(k)}  be  two  sequences.  If  there  exists  a sequence 
{^(k)}  with  8(k)  0 as  k -►  » such  that  y(k)  - 3(k)x(k)  then  we  denote 

y(k)  - o[x(k)] 

Definition  2;  Let  {x(k)}  and  {y(k)}  be  two  sequences.  If  there  exists  a positive 
constant  M such  that 

|y(k) I J M|x(k) I for  all  k e N 

We  denote  It  by 

y(k)  - 0[x(k)] 

Definition  3;  If  {x(k)}  and  {y(k)}  are  two  sequences  such  that  x(k)  - 0[y(k)] 
and  y(k)  > Olx(k) ] , we  say  that  the  two  sequences  are  equivalent  and  denote  It  by 

x(k)  -v  y(k). 

The  following  definition  Is  found  to  be  useful  while  describing  two  sequences 
which  evolve  at  the  same  rate  but  are  not  equivalent. 

Definition  4;  Two  sequences  {x(k)}  and  {y(k)}  are  said  to  grow  at  the  same  rate  If 

sup|x(v)  I “V  sup|y(v)  I 
k^v  k^v 

It  follows  that  two  sequences  which  are  equivalent  grow  at  the  same  rate  but 
not  vice  versa. 

Definition  5t  Sequence  (x(k)}  Is  said  to  grow  faster  than  sequence  {y(k)}  If 

sup|y(v)|  - o[sup|x(v)|] 
k^v  k^v 

Let  V^(s)  be  a rational  transfer  fiinctlon  of  a linear  tlae-lnvarlant  discrete 
system  with  all  Its  poles  and  seros  within  the  unit  circle  and  Input  and  output 
x(k)  and  y(k)  respectively.  Let  h(k)  be  the  impulse  response  of  W^(e).  Since 


f 


r 


I 
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i 


h c It  follows  that 


ly(k) I ^ c 8up|x(v) I + c, 
^ k>v  ^ 


(A.l) 


where  and  are  positive  constants.  Hence 


8up|y(v)|  - 0t8uplx(v)|] 
kjv  kiv 


,-1. 


Considering  y(k)  as  the  input  and  x(k)  as  the  output  we  also  have  (since  (z) 
is  stable) 


|x(k) I ^ c'  sup  |y(v) ( + c' 
A 1^—  ^ 


(A.  2) 


k+rj^iv 


where  c^  and  c^  are  positive  constants  and  r^  is  the  relative  degree  (i.e.  number  of 


poles  - number  of  zeros)  of  U (z).  If  the  rate  at  which  the  sequence  y(k)  can  grow 


is  bounded  [e.g.  any  linear  system  with  bounded  coefficients]  it  follows  from  A. 2 
that 

(A.3) 


(x(k)(  ( c"  8up|y(v)|  + c' 
k^v 

*^1 

"1  "2^ 


cV  is  a positive  constant  and  c"  c'|xj  where  X.  denotes  the  maximum  rate  at  which 


y(k)  can  grow. 

In  view  of  (A.l)  and  (A.3) 


sup|x(v)|  “V  sup|y(v)| 
k^v  k^v 


(A.  4) 


and  {x(k)}  and  {y(k)}  grow  at  the  sasw  rate. 


Let  H (z)I  be  a transfer  matrix  with  input  vector  w(k)  and  output  vector  ((k) 
L 


where  H.  (z)  is  a minimum  phase  transfer  function  and  I is  the  unit  matrix.  If  w(k) 

is 


and  ((k)  are  repreeenced  as 


2n 


w(k)  - E Wj^(k)£. 
i-1  ^ 

2n 

C(k)  - E C4(W)l. 
i-1  ^ 


(A.5) 


t ! 


I ■ 
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where  are  Independent  constant  vectors  it  follows  that  and 

froB  (A. 4)  that 

aupU  (v)l  'V.  8up|w.(v)  1 (A. 6) 

k^v  ^ kjv 

1 - 1,2, ...2n 
2n 

Definition  6;  Let  c(k)  - I c.(k)l.  where  I are  Independent  constant  vectors. 

1-1 

The  components  C^(k)  (1  - 1,2 are  said  to  be  dominant  if  Cj^(k)  grows  at 

the  same  rate  as  ||c(k)||  or 


®'iplc.(v)|  'v«  supi|c(v)||  1 - 1,2, ...m  (A. 7) 

IttV  ^ kjv 

By  definitions  (5)  and  (6)  the  components  Cj^(v)  1 - m^^+l 2n  would  grow  at  a 

slower  rate  than  |c(k)||  and 


supl4.(v)|  - o[8up  ||  C(v)||  1 1 - m^+l,...,2n  (A. 8) 

kiv  ^ kiv 


Lou:  If  w(k)  and  C(k)  ate  the  vector  Input  and  output  respectively  of  a minimum 


phase  transfer  matrix  U (z)I  such  that  (A. 4)  holds 

L 

and  w^(k)  1 - l,...»mj^  are 

dominant  then 

sup|w.(k) 1 

-V,  auplU(v)  II 

1.  “ l.p29e«e 

k^v 

k;v 

(A.  9) 

eup|w  (k) 1 

- o(oup  (|c(v)  II  1 

i ■ m^+1, . . .2n 

kiv 

kiv 

Proof : Follow 

directly  from  (A.5) ,(A.6) ,(A.7)  and  (A. 8). 

Further  since 

sup  ||w(v)|| 

sup  II  i:(v)|| 

kiv 

k^v 

we  have 

eup|w  (v) 1 

aupIlwCv)!! 

1.  *■  Lp  2p  • • eft^ 

k>v 

k)v 

end 

eup|w.(v) 1 
kbv  ^ 

- o[sup  |w(v)  II  ] 

1 t • • • 1 2x1 

(A.  10) 

so  that  w^(k) 

i ■ l,2,...,m.  are  dominant. 

2n  ^ 

Definition  7t 

If  C(k)  - I C.(k)t.  with  C.(k)  1 

i-1  ^ ^ ^ 

• lp2p...»a^  doalttibntp  th«n 

^^(k) 

are  aald  to  be  asymptotically  independent  if 
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S c.C.(k)  - o[8up  ||j;(v)||  ] 
i-1  ^ ^ tev 


(A. 11) 


Implies  c,  ■ c*  ■ ...  ■ c ■ 0 
X X 

without  loss  of  generality  we  can  always  assume  that  a number  m^, 

1 J m,  ^ 2n  exists  such  that  c,(k),...,C  (k)  are  asymptotically  Independent. 

1 ^ ®'l 

If  c^(k)  (1  ■ l,2,...,m^)  are  asymptotically  Independent  and 

“l 

E c (k)C.(k)  - o[8up  ||c(v)||  1 
1-1  ^ ^ kiv 

and  11m  c^(k)  exists  for  all  1 - l,2,...,nL  then  It  follows  that  11m  c (k)  - 9 for 
^ ^ kx»  ^ 

all  1 ■ l,2,...»m^. 

From  equations  (10), (11)  and  (12) 

♦'^(k)c(k)  - aJ(k)r;(k)ej^(k)  -0 
♦ (k)  -*•  ♦* 


C(k)ej^(k)  -►  0 


as  k -»■  «>. 


Hence  ac^(k)rc(k)ej^(k)  - oI|jc(k)||  ] and  It  follows  that 


U’^(k)c(k)l-  o[||c(k)||  ] - ot8up||i:(v)||  ] 

k^v 


If  Cj^(k)  (1  ■ l,2,...,mj^)  are  asymptotically  Independent 


♦ - 0 


1 — 1,2,... ,m- 


(A. 12) 


By  (A. 10)  and  (A. 12)  It  follows  that 


sup|^^(k)w(k)  I - o[8up  ||w(v)||  ] 
k^v  k^v 


(A. 13) 


